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Phase transition of charged black holes in massive gravity through new methods 
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Motivated by providing preliminary steps to understand the conception of quantum gravity, in 
this paper, we study the phase structure of a semiclassical gravitational system. We investigate the 
stability conditions and phase transition of charged black holes in massive gravity via canonical and 
geometrical thermodynamic approaches. We point out the effects of massive parameter on stability 
conditions of these black holes and show how massive coefficients affect the phase transition points 
of these black holes. We also study the effects of boundary topology on thermodynamical behavior 
of the system. In addition, we give some arguments regarding the role of higher dimensions and 
highlight the effect of the electric charge in thermodynamical behavior. Then, we extend our study to 
geometrical thermodynamic approach and show that it can be a successful method for studying the 
black hole phase transition. At last, by employing the relation between thermodynamical pressure 
and cosmological constant, critical behavior of the system and the effects of different parameters on 
critical values are investigated. 


I. INTRODUCTION 


Cosmological observations show that, 95% of the Universe is made of dark energy and dark matter. In order to 
interpret dark sector of the Universe, some various candidates have been proposed in literature, such as modified 
gravities like Lovelock gravity M, brane world cosmology scalar-tensor theories [I® , and the so-called 

F(R ) gravity theories [U® (for a review regarding thermodynamical aspects of modified theories, see Ref. i3)- 
One of the interesting aspects of these theories is graviton. Graviton in the Einstein gravity is massless A 

natural question is whether one can build a self-consistent gravity theory if graviton is massive (4ll-l43j. To do so, one 
should modify the Einstein theory of gravity to include massive gravitons ill I hsl . From the theoretical perspective, 
the shear difficulty of constructing a consistent theory of massive gravity makes the problem more interesting. 

Fierz and Pauli developed the ghost-free theory of non-interacting massive gravitons which had a condition of being 
linearized over Minkowski space [49|,|5C|. Then, Boulware and Deser showed that a generic extension of the Fierz-Pauli 
theory to curved backgrounds will contain ghost instabilities [51,]. In Refs. [46l. |48[ [52 - l55[ . a potentially ghost-free non¬ 
linear massive gravity actions has been constructed. Furthermore, one of the classes of charged black hole solutions in 
nonlinear massive gravity was studied in Ref. [5f|. In addition, the effective field theory of the modified gravity with 
massive modes has been studied in Ref. [57j ■ The phenomenological aspects of massive gravities with Lorentz-violation 
and invariant mass terms were investigated 58j- Furthermore, massive, massless and ghost modes of gravitational 
waves in higher order gravities were studied and the possibility of the detection and additional polarization modes 
of a stochastic gravitational wave were discussed in Ref. j59|. In a series of studies, the black hole solutions in 
different massive gravities have been investigated and it was shown that some of these black holes may enjoy 

anti-evaporation in their thermodynamical properties [7ll. f72l]. The massive gravity has also been employed to conduct 
some studies in the context of cosmology (for an incomplete list of these studies, we refer the reader to Refs. [73i48lj l. 
Of the other applications of massive gravity, one can point out the large applications of this gravity in the anti-de 
Sitter/conformal field theory (AdS/CFT) community. For example, the lrolographical aspects of the massive gravity 
coupled with nonlinear electrodynamic models have been investigated in Ref. [82| . Furthermore, the extension of 
massive gravities with additional couplings and higher derivative terms have been addressed and their holographical 
properties are extracted 83 h86|. 

Recently, David Vegh |87| found a nontrivial black hole solution with a Ricci flat horizon in four dimensional massive 
gravity with a negative cosmological constant. He studied an assumption that the graviton may have lattice like 
behavior in the holographic conductor model: the conductivity generally exhibits a Drude peak which approaches to a 
delta function in the massless gravity limit. Application of massive gravity in the context of AdS/CFT correspondence 
via stability conditions and metal-insulator transition has been studied in [88l - l9i| . In addition, thermodynamical 
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behavior and P—V criticality of the Vegh’s massive gravity have been, recently, investigated I92H961]. Some holographic 
effects of graviton mass have been investigated in Ref. [93, j98( . The main goal of this paper is to obtain Vegh’s black 
hole solutions and study corresponding thermodynamical properties and their phase structure. 

Thermodynamical aspects of the black holes and their critical behavior have been one of the greatest interests 
for many authors. Recently, it was proposed to change the role of some constants to the cases of (dynamical or 
thermodynamical) variables f99t-Tl07i|. Among the thermodynamical properties of the black holes, thermal stability 
and phase transition have been investigated in literature. Studying thermal stability and phase transition points can 
be done in the context of canonical ensemble by analyzing the heat capacity. The sign of heat capacity represents 
stability/instability of the black holes, while its divergencies and roots are denoted as second order phase transition 
and bound points, respectively. 

Another approach toward studying critical behavior of the black holes is through geometrical method. In this case, 
by employing a thermodynamical potential and corresponding extensive parameters, one can construct a phase space 
containing information regarding phase transition of the black holes. In other words, geometrical thermodynamics 
tries to introduce thermodynamical properties of the systems through the use of the Riemannian geometry. The 
basic motivation is to give an independent picture regarding thermodynamical aspects of systems. In addition, 
geometrical thermodynamics gives information regarding bound points, phase transitions, their types and stability 
conditions. Furthermore, it contains information regarding molecular interaction around phase transitions for the 
mundane systems. In other words, by studying the sign of thermodynamical Ricci scalar (TRS) around phase transition 
points, one can extract information whether interaction is repulsive or attractive. The re are sev eral approaches in the 
context of geometrical thermodynamics (GTs); the well-known metrics of the Weinhold Il08lll09| . Ruppeiner [01 on 
and Quevedo metric which were proposed in different structures. Recently it was pointed out that ther e may be some 
mismatches between second order phase transitions and divergencies of TRS in the mentioned methods [ll2Mll5 | . and 
an alternative thermodynamical metric was introduced to solve this problem eimh. The generalization__of this 
metric in the context of extended phase space (in which constants considered as variables) was done in Ref. |l20i | . It 
was shown that thermodynami cal b ehavior of the system in the context of extended phase space, heat capacity and 
GTs lead to same consequences |l2dj |. The importance of geometrical methods comes from the fact that these methods 
are based on the grand canonical ensemble foundation. Although we connect this method to the heat capacity, it is 
worthwhile to mention that this is only for the sake of comparison and check the validity of the results. 

Recent progresses in the Ad S/CF T correspondence, suggest that cosmological constant should be considered as a 
thermodynamical variable jl2lMl24 |. In addition, in black hole thermodynamics, it wa s sho wn that such consideration 
could lead to remove specific thermodynamical problems such as ensemble dependency jl25| . One of the corresponding 
thermodynamical variable for cosmological constant is pressure. It was shown that by using this analogy, one can 
introduce novel properties in black hole thermodynamics such as Van der Waals like behavior, triple point and 
reentrant of phase transition 12d4l35| |. Therefore, by such consideration, it is possible to study critical behavior of 


the black holes in extended phase spa ce. One of the new methods for extracting critical points of the black holes was 
proposed and employed in Ref. [120 ], In this method, by using proportionality between cosmological constant and 
pressure in denominator of the heat capacity, a relation is obtainable for pressure which is independent of equation 
of state. In this paper, we will employ this method to extract critical points of the black hole system. 

The outline of the paper is as follow; in section II, we review the charged massive black holes and their thermody¬ 
namical quantities. In section III, we introduce approaches for studying phase transition of these black holes in the 
context of heat capacity and geometrical thermodynamics, and study thermal stability of these black holes. Then, 
we investigate the existence of the phase transitions in the context of two mentioned approaches. The last section is 
devoted to closing remarks. 


II. BLACK HOLE SOLUTIONS IN MASSIVE GRAVITY 


The (n + 2)-dimensional action of massive gravity with negative cosmological constant can be written as 

n+2 
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where 1Z is the scalar curvature, T is t he M axwell invariant, / is a fixed rank-2 symmetric tensor and m 2 is the 
positive massive parameter (see Ref. j|87j, Il36j for more details regarding the relation of positive m 2 with the wall of 
stability interpretation). In addition, Ci s are constants and lAi s are symmetric polynomials of the eigenvalues of the 
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(n + 2) x (n + 2) matrix K.% = y/g^fa. 

Ui = [AC], 

= [/C] 2 - [/c 2 ], 

U 3 = [/C] 3 - 3 [AC] [/C 2 ] + 2 [AC 3 ] , 

U 4 = [/C] 4 - 6 [/C 2 ] [AC] 2 + 8 [AC 3 ] [AC] + 3 [AC 2 ] 2 - 6 [AC 4 ] , 


U 3 = [AC] 5 - 10 [AC] 3 [AC 2 ] + 20 [AC] 2 [AC 3 ] - 20 [AC 2 ] [AC 3 ] + 15 [AC] [AC 2 ] 2 - 30 [AC] [AC 4 ] + 24 [AC 5 ] , 


l 2 rir3 


( 2 ) 


The square root in AC means an d [AC] = AC]]. Variation of the action JT} with respect to the 

metric tensor and the Faraday tensor F^ v , leads to 


G, lv - 


i(n + 1) 1 

2Z 2 9 ^ = 2 




m 2 x^, 


(3) 


V, t F^ = 0, 


(4) 


where G is the Einstein tensor and 

x, jv = -y {Uigp, - AC^) - y (i U 2 g W - 2U[ AC,,,, + 2AC 2 J - y (U 3 g^ - 3ZY 2 AC^ + 

6WxAC 2 , - 6/C 3 .„) - y (U 4 g^ - \U 3 KL [lv + 12W2AC 2 , - 242V,ACf_ + 24AC 4 J + ... . (5) 

We take into account the metric of (n + 2)-dimensional spacetime with the following form 

ds 2 = —f(r)dt 2 + f^^dr 2 + r 2 hijdxidxj, i, j = 1,2,3, ...,n (6) 

where hijdxidxj is the line element for an n —dimensional Einstein space with constant curvature n{n + 1 )k and 
volume V n . We should note that the constant k indicates that the boundary of t = constant and r = constant can 
be a positive (elliptic), zero (flat) or negative (hyperbolic) constant curvature hypersurface. 

A generalized version of / M „ was proposed in Ref. I92HD5I] with the following form 


fiiv = diag{ 0 , 0 , clh tj ), 

where by employing it, the metric function f(r) will be [92i Tf)^] 

r 2 

Z 2 


f(r) = k +--^ + _ Q - _+ CoClw2 r - 

12 r n i 2 n(n — l)r 2 (”- 1 ) n 


2 2 
c 0 c 2 m 


(n — l)cgC 3 m 2 r 1 + (n — l)(n — 2 )cQC±m 2 r 2 + ... 


(7) 


( 8 ) 


in which too and q are integration constants which are related to the total mass and electric charge of the black hole, 
respectively. As one can see, the massive terms could be up to n + 2 number. Technically, it is not easy to study the 
total behavior of the solut ions with this number of terms. Therefore, for the simplicity, we restrict our study to Ui 
up to the fourth term, U 4 [l53l |. 

In order to study the effects of massive gravity, one can investigate the metric function. Regarding various terms of 
/(r), it is worthwhile to mention that fourth term ( q 2 term) is dominant near the origin and therefore one can conclude 
that the singularity is timelike. In addition, for large distance, second term (r 2 /l 2 term) is dominant which confirms 
that the solutions are asymptotically adS. It is evident that the behavior of the metric function is highly sensitive 
to massive parameters and contribution of the massive gravity (see Figs. [2] and [3]). For specific values of different 
parameters, the total behavior of the metric function will be similar to the metric function of Einstein-Maxwell 
gravity which is containing two roots, one extreme root or without any root (see Fig. [T]). Whereas by considering 
other choices of free parameters, the following configurations for the roots of the metric function may happen. For 
example, two extreme roots (dashed line in Fig. [2]), three roots in which the smaller one is extreme (dotted line in 
Fig. [2]), four roots (bold-dotted line in Fig. [2]), three roots where the bigger one is extreme (continues line in Fig. 
m and three roots which the middle one is extreme (bold continues line in Fig. [2]). As one can see, these different 
behaviors are originated from contributions of the massive gravity. The existence of multi-horizon for massive gravity 
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FIG. 1: f(r) versus r for l = 1, q — 1, Co = 0.1, ci = —0.3, C 2 = 2, C 3 = 3, C 4 = 2, fc = 1 and n = 4. 

Left diagram for mo = 0.5, m = 0.3 (dashed line), m = 2.38 (continues line) and m = 3 (dotted line). 

Right diagram for m = 2, mo = 0.45 (dashed line), mo = 0.49 (continues line) and mo = 0.55 (dotted line). 



FIG. 2: f(r) versus r for l = 1, q = 1, Co = 0.1, ci = —2.3, C 2 = 10, k = 1, and n = 2. 

Diagrams for mo = 1.648, m = 5.50 (dashed line), mo = 1.668, m = 5.60 (dotted line), mo = 1.70, m = 5.65 (bold-dotted line), 
mo = 1.68, m = 5.462 (continues line) and mo = 1.74, m = 5.65 (bold line). 


has been reported in Refs. 0,113 as well. It was shown that the existence of such property with specific quantum 
corrections leads to a phenomena which is opposite of the black holes evaporation. In other words, there might exist 
an anti-evaporation behavior for the black holes in these cases. Here, we see that similar property is observed for 
these black holes. Therefore, it is possible that anti-evaporation property exists for these black holes in the presence 
of massive gravity. I t is worth while to mention that anti-evaporation has been reported for specific black holes in 
F(R) gravity as well fTMITifl . 

Considering the fact that we are interested in studying GTs and heat capacity of the solutions, we give a brief 
















5 




r r 




r r 

FIG. 3: f(r) versus r for l =s1, q = 1, co — 0.1, C 4 = 1, k = 1, mo = 0.6, m = 5.5 and n = 4. 

Left-up panel: for C 2 = C 3 = 1, ci = —1.7 (continuous line), ci = —1.44 (dotted line), ci = —1.3 (dashed line) and c\ = —1.23 
(dashed-dotted line). 

Right-up panel: for C 2 = C 3 = 1, ci = —0.5 (continuous line), ci = 0 (dotted line), ci = 1.46 (dashed line) and ci = 2 
(dashed-dotted line). 

left-down panel: for ci = —1.7, C 3 = 1, C 2 = 1 (continuous line), C 2 = 1.78 (dotted line), C 2 = 4 (dashed line) and C 2 = 4.85 
(dashed-dotted line). 

right-down panel: for ci = —1.7, C 2 = 1, C 3 = 3.8 (continuous line), C 3 = 6 (dotted line), C 3 = 8.5 (dashed line) and C 3 = 12 
(dashed-dotted line). 


review about the conserved charges and thermodynamic quantities in which they have been obtained in Ref. [92| ■ 
Regarding the obtained solutions, the action © diverges, as is the Hamiltonian and other associated conserved 
quantities. One of the systematic methods f or re moving this divergency is through the use of the counterterm method 
inspired by AdS/CFT correspondence Il4ll4l44| . It was also shown that by using the Hamiltonian approach, one can 
find the mass M of the black hole solutions as ( see Ref. |92j for more details) 

M = ( 9 ) 

In addition, the electric charge of the black hole, Q, can be obtained by using the Gauss’s law and calculating the 
electromagnetic flux at infinity. Straightforward calculations lead to the following result 


Q = 2^“- 


( 10 ) 


The temperature may be obtained through the use of regularity of the solutions at the event horizon (the largest 
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root of metric function which is denoted by r+). It is also calculated through the use of the surface gravity definition 
as 


1 


1 


27r 

where % = dt is the Killing vector. It is easy to show that 

1 


T= ^-\—z{VnXv)('W*X v ), 


T = 


47rr_i 


r 2 o 2 

(n- l)fc + (n+ , .. +cic 0 m 2 r + (n- 1 )c 2 CqTO 2 + 

1 2nr2 J 


(n — l)(n — 2 )c 3 CqTO 2 (n — l)(n — 2 )(n — 3)c4CqTO 2 


r+ 


( 11 ) 


where r + satisfy f(r = r+) =0. 

At last, we should o btain the entropy of the black hole solutions. The entropy of the black holes satisfies the 
so-called area law [l45l - ll5(ij in the context of Einstein gravity. According to this law, one finds that the black hole 
entropy is equal to one-quarter of the horizon area, i.e., 


S = 


2nV n 

K 2 


'+■ 


( 12 ) 


III. GEOMETRICAL STUDY OF THE PHASE TRANSITION AND THERMAL STABILITY 

In this section, we study the local stability and second order phase transition of the solutions in the canonical 
ensemble and compare our results with a new GTs approach for the phase transition and bound points. 

In the canonical ensemble, the positivity of the heat capacity is sufficient to ensure thermal stability. One can 
calculate the heat capacity as 


r _ r (dS\ ( os)q , . 

Cq - t [ Trr ~ (a 2 M\ ■ ( 13 ) 

\ a± ' Q \-dsr) Q 

Considering the fact that changing in the sign of heat capacity is representing the phase transition between unsta¬ 
ble/stable states, we regard divergence points of the heat capacity as second order phase transition points. So, the 
second order phase transition and bound points of the black holes with regular T in the context of heat capacity are 
indicated with following relations 


T = (Off) q = 0 bound point 


( d 2 M \ 

vw) 


Q 


= 0 phase transition point 


(14) 


where we called the roots and divergence points of the heat capacity as bound points and (second order) phase 
transition points, respectively. It is notable that T = 0 indicates a bound point between nonphysical (T < 0) and 
physical (T > 0) region. 

On the other hand, GTs is another way for investigating the phase transition in the context of black hole ther¬ 
modynamics. In this approach, several metrics have been introduced (Ruppeiner, Weinhold and Quevedo metrics) 
in which by using of these metrics, one can investigate the phase transition of black h oles. It was previously shown 
that these metrics encounter with some problems for specific types of black holes [l!6l - lll9l | . R ecent ly, a new metric 


(HPEM metric) was proposed in order to solve the problems that other metrics may confront [nMii. The roots 


of denominator of the Ricci scalar of HPEM metric only contains the roots of numerator and denominator of the 
heat capacity. In other words, divergence points of the Ricci scalar of HPEM metric only coincide with the roots and 
divergence points of the heat capacity. The HPEM metric has the following form 


ds 2 

ab HPEM 


SMg 


(TTn d 2 M S \ 

V i=2 dxi J 


-MgsdS 2 


i —2 


/ d 2 M 


dXi 


(15) 
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where Ms = dM/dS, Mss = d 2 M/dS 2 and \i (Xi 7^ >5) are extensive parameters. In what follows, we will study 
the stability and second order phase transition of the charged massive black holes in the context of heat capacity and 
GTs. 

First of all, one should take this fact into consideration that the sign of temperature is putting a restriction on 
systems to be physical or non-physical. As we see later, there will be a critical horizon radius, r+ c , in which for 
r+ < r+ c the temperature of the system is negative and in this region solutions are non-physical. 

It is a matter of calculation to show that for the root of the heat capacity, one can find following relation with 
respect to the massive parameter 


where 


— I 1 [A 

m c =m \ Cq=q = -^-, 


(16) 


A = 2 \q 2 l 2 — 2 n (n + 1) r+ n ] r+ — 4n/ 2 /cr+™ 4 , (17) 

B = nc 0 r 2 ” -4 [c 0 c 2 (n - 1) r+ + cir+ + (n - 1) (n - 2) eg {r+c 3 + c 0 c 4 (n - 3)}] . (18) 

Regarding positive c* leads to positive B and therefore, we should restrict the parameters to obtain positive A. The 
limitation reduces to the following inequality 

q 2 l 2 - 2 knl 2 r 2 +~ 5 -2 n(n + 1) > 0, 

where for the flat boundary (k = 0), one finds 


J q 2 l 2 

+ \2n(n+l)J ' 

It is worthwhile to mention that in the absence of massive parameter (m = 0), the obtained solutions will reduce 
to Reissner-Nordstrom black holes and these black holes only enjoy bound point. Here, in the presence of massive 
gravity (m ^ 0), there is also only one root for heat capacity, which is a function of the massive parameter. It is 
notable that massive coefficients (cj’s), are only presented in denominator of the obtained critical value. Therefore, 
in case of Cj > 1, the effects of these coefficients on the critical value of massive parameter, are of decreasing ones. 
Whereas for the case of 0 < Cj < 1, their effects are in favor of increasing the value of the critical massive parameter. 
In other words, they increase the value of m c - 

Next, as for the divergence points of the heat capacity, hence second order phase transitions, we find following 
critical relation for the massive parameter 


in which 


1 

fc 


’"• sra A-“ = 2i^\ 

' D' 

(19) 

C = 6 [(2n — 1) q 2 l 2 + 2n (n + 1) r 2n ] r+ — 

12kl 2 n(n-l)r 2 ™~ 2 , 

(20) 

D = 3n (n — 1) r 2n 4 [r+c 2 + c 0 (re — 2) {2r+c 3 + 3c 0 C4 (n — 3)}] . 

(21) 


It is evident that in this case, similar to the case of heat capacity’s root, due to structure of the square root 
functions, there will be some restrictions. Interestingly, in this case for k = 0 and k = — 1 and due to our interest 
only in positive values of massive coefficients, C and D are always positive and there is no restriction. Whereas for 
the case of spherical symmetric, there is only one restriction which is in the following form 

(2 n - 1) q 2 l 2 + 2 n (n + 1) r 2n - 2 l 2 n (n - 1) r ^ n ~ 3 > 0, (22) 


in which due to its complexity it was not possible to find an analytical restriction with respect to horizon radius. 
Another interesting property of the obtained equation for the divergence points of the heat capacity is the independency 
of the critical massive parameter of the c 4 coefficient. In other words, critical value of the massive parameter in case 
of divergencies of the heat capacity is independent of the variation of c 4 . Therefore, there is no contribution of the c 4 
to critical value of the massive parameter. It is worthwhile to mention that, in this case similar to root of the heat 




FIG. 4: rrii (continues line) versus r+ for l = q = 1, Co = C 2 = C 3 = C 4 = 2, n = 4, and k = 1. mi = 0.04 (dotted line) mi = 0.4 
(dashed line) 
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FIG. 5: 1Z (continuous line), Cq (dashed line) and T (dotted line) versus r+ for l = q = 1, co = ci = C 2 = C 3 = C 4 = 2, n — 4, 
and k = 1 ; m = 0 (up - left), m = 0.01 (up - right) and m = 1 (down panels with different scales). 
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FIG. 6: TZ (continuous line), Cq (dashed line) and T (dotted line) versus r+ for l = q = 1, ci = C 2 = C 3 = a ## 2, n s= 4, 
m = 0.1 and k = 1; Co = 0 (up - left), Co = 1 (up - right) and Co = 2 (down panels with different scales). 


capacity, the presences of the massive coefficients are only observed in denominator of mt. Therefore, similar effects 
for variation of the massive coefficients in case of root of heat capacity will be observed in this case too. 

In addition, it must be pointed out that the obtained value of massive parameter in Eq. (1191) is not an injective 
function. Plotting nii with respect to horizon radius, one finds that there is a value of rrii, say rrii c in which for 
rrii < nii c , there will be no critical horizon radius available. In other words, in this region there is no divergence point 
for the heat capacity. In case of rrii = rrii c , there will be only one divergence point for heat capacity and finally for the 
case of rrii > rrii c , there will be two horizon radii for any critical value of the massive parameter. This property of the 
obtained equation will be shown in the plotted graphs for the heat capacity. Here, in order to elaborate mentioned 
behavior for Eq. m, we have plotted Fig. [4j 

Next, by employing the Eqs. (mu, m and ED, we plot Figs. \m to study thermodynamic and geometrical 
thermodynamic behavior of these black holes. 

It is evident that for the small values of massive parameter, there is no second order phase transition (Fig. [5] 
up). Whereas for large values of this parameter, there is a bound point and a second order phase transition (Fig. [5] 
down). It is worthwhile to mention that this behavior could be observed by studying Fig. [U Also, according to the 
plotted graph for to,; (Fig. Hj), the smaller (larger) divergence point of the heat capacity is a decreasing (an increasing) 
function of m. This behavior is also evident from plotted graphs for variation of m (Fig. [5]). 

Besides, due to coupling of different powers of cq with massive parameter, we have studied the effects of variation 
of this coefficient. Interestingly, similar to behavior of the massive parameter, for small values of Co, only a bound 
point is observed (Fig. [G]up). It is seen that only for sufficiently large values of this coefficient, two second order 
phase transitions and a bound point are obtained (Fig. [G]down). 

Our next effects of interest are due to variation of topology of the solutions. If we denote the indexes F, H and S 
as representing flat (k = 0 ), hyperbolic (k = — 1 ) and spherical horizons respectively, and also root of heat capacity 
with ro and two divergence points of the heat capacity with r^ivi (smaller one) and roiv 2 (larger one), we can find 
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FIG. 7: 1Z (continuous line), Cq (dashed line) and T (dotted line) versus r+ for l = q = 1, Co = ci = ci = cz = ca = 2, n = 4, 
and m = 0.1; for different scale: k = 1 (left panels up and down), k = —1 (middle panels up and down) and k = 0 (right panels 
up and down). 




FIG. 8 : Cq versus r+ for l = q = 1, Co = Ci = C 2 = C 3 = C 4 = 2, k = 1, and m = 0.1. 

for different scales: n = 4 (continues line), n = 5 (dashed line), n = 6 (dotted line) and n = 7 (bold continues line). 
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FIG. 9: m,i (continues line) versus r+ for l = 1, Co = ci = C 2 ~ C 3 = C 4 = 2, k = l, n = 4 (continues line), n = 5 (dashed line), 
n = 6 (dashed-dotted line) and n = 7 (dotted line); q = 1 (left panel) and q = 0 (right panel) 


following results. For the root of the heat capacity, we have ro_s < tq< r^-H- hr other words, the highest 
value for the root of heat capacity belongs to the solutions with hyperbolic topology. In addition, one can find that 
fDivi-s < roivi-H < rDivi-F which indicates that the highest value of smaller divergence point belongs to flat 
horizon (Fig. [7] up). Finally, for the larger divergence point, one can find that rr>ivi-H < rDivi-F < rDivi-S which 
means that in this case the highest value of larger divergence point belongs to spherical symmetric black holes (Fig. 
[7] down). 

Finally, we are in a position to study the effects of dimensions on thermodynamical behavior of the system. In 
plotted graph, we only present the heat capacity diagrams for various dimensions. Evidently the root of the heat 
capacity is only an increasing function of dimensions (Fig. [5] left). As for divergence point of the heat capacity, the 
larger one is also an increasing function of dimensions (Fig. [5] right). But regarding the smaller divergence point, 
interestingly, an anomaly is observed in case of n = 5. In other words, for increasing value of dimensions, the smaller 
divergence point increases except for the case of n = 5 which its smaller divergence point is smaller than the case 
of n = 4. In order to find the reason for this behavior, once again we plot Fig. [9] (left) for different dimensions. 
Evidently, for specific small values of m , the behavior of the system is with order, in which increment in dimensions 
leads to increase in value of the smaller divergence point of the heat capacity. On the other hand, for special values 
of to the smaller divergence points of different dimensions will be the same. In other words, there are some values of 
to in which smaller divergence point of the heat capacity will be the same for n = 4 and n = 5 or other dimensions. 
This behavior of the system is solely because of the structure of to,, and it is the characteristic behavior of the system 
with this configuration (for this type of massive gravity). It is worthwhile to mention that increasing the value of 
dimensions will decrease the region in which there is no second order phase transition. In other words, the minimum 
of the mi is a decreasing function of dimensions. Interestingly, for the case of vanishing electric charge (Fig. fright), 
the mentioned behavior will not be seen. In this case, there is only one second order phase transition for these black 
holes which specifies the role of the electric charge in thermodynamical behavior of the system. 

Now, we focus to study thermal stability of the solutions. The stability of the solutions, hence, thermodynamical 
behavior of these black holes are presented by the bound point and number of the second order phase transitions. For 
example, in case of only one bound point (up panels of Figs. [5] and [7]) , there is a region in which both temperature 
and heat capacity are negative. Therefore, this region is presenting the non-physical unstable state for black holes. 
After passing the bound point, the system acquires a stable physical state. For the cases of two second order phase 
transitions and one bound point (down panels of Figs. [G]and[7]), there are several changes in the sign of heat capacity. 
By using naming that was used before, one can find that for the case of r + < ro, system is in nonphysical and 
unstable phase. Whereas in r+ = r 0 system goes under a transition from nonphysical to physical state and for 
ro < r_|_ < roiv 1 , system is in thermally stable state with positive temperature. As for the roiv 1 < r+ < rm v 2 , system 
is in unstable state and for rm v 2 < r+, black holes are thermally stable with positive temperature. Therefore, due to 
thermodynamical concept in which system desires to find stable state, there are two second order phase transitions. 
For the case of r + = roiv 1 + <5, hr which S is so small, the system goes under a second order phase transition from 
larger unstable black hole to smaller stable black hole. In other words, in this phase transition the size of the black 
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hole decreases. On the other hand, for the case of r + = rmvi — S, black hole in this configuration goes under another 
second order phase transition. In this one, black hole acquires larger horizon radius, hence its size increases. 

Next, we study the geometrical thermodynamic behavior of the system. Figs. [5]to[8]show that considered metric in 
the context of GTs, provided a successful machinery. In other words, bound point and second order phase transitions 
for the heat capacity and the divergence points of Ricci scalar of the considered metric coincide with each other. 
Therefore, the mentioned metric has divergencies in place of bound point as well as second order phase transition 
points. On the other hand, the behavior of the Ricci scalar near bound point and second order phase transitions is 
not the same. For the root of the heat capacity (bound point) the sign of TRS is different for before and after that 
point. Whereas, in case of the second order phase transition, the sign of TRS is fixed. This change in sign enables 
one to distinguish bound point from second order phase transition. Therefore, it is possible to find bound point and 
phase transition by applying GTs, without the use of heat capacity diagrams. 


IV. CRITICAL BEHAVIOR IN EXTENDED PHASE SPACE 

A. critical behavior in extended phase space through usual method 

The van der Waals like behavior of the black holes could be studied through the use of the analogy between cos¬ 
mological constant and thermodynamical pressure. The proportionality between cosmological constant and pressure 
is given by 


P = - 


A 

87T 


(23) 


By replacing the cosmological constant with thermodynamical pressure in temperature, one can obtain the equation 
of state. Using this equation of state and inflection point property, one can extract critical points and study van der 
Waals like behavior. For black hole solutions which considered in this paper, the critical behavior of the system 
through usual method has been investigated in Ref. [933. Here, we would like to investigate the critical properties of 
the black holes through a new method. 


B. critical behavior in extended phase space through new method 

In this section, we will study the critical behavior of these black holes through the use of the prop ortio nality between 
pressure and cosmological constant. The method which we employ to do so, was introduced in Ref. [l2(| . By replacing 
the cosmological constant with its corresponding thermodynamical pressure in the heat capacity, one is able to obtain 
a relation for the pressure by solving the denominator of heat capacity. The maximum of P(r+) is where the second 
order phase transition takes place. In this relation, P does not explicitly depend on the temperature, and therefore, 
this relation differs from equation of state, but it contains all the information regarding phase structure of the black 
holes. For the pressures smaller than critical pressure, two horizon radii for each pressure exist which indicate three 
phases. While for pressures larger than critical pressure, no phase transition is observable. 

Using the denominator of heat capacity (fl3l) with (l23l) . one can find following relation for the pressure 

3n(n — l) 2 (n — 2)(n — 3)cgm 2 C4 n(n — l) 2 (n — 2)cgm 2 C3 n(n — 1) 2 CoTO 2 c 2 
167rr + (l + n) 87rrj.(l + n) 167rr+(l+n) 

(n- l)(2re- l)q 2 kn(n - l) 2 
327rr^_ n (l + n) 16nr+(l + n)' 

In order to find the critical horizon radius, one should calculate first order derivation of the pressure with respect 
to horizon radius and solve it for r+. Calculations show that it is not possible to obtain critical horizon radius 
analytically. Therefore, we employ numerical method. The results are presented in various diagrams (Figs. I10H12I) . 
In addition, we have plotted coexistence curve which characterizes two different phases of small/larger black holes 
with same temperature and pressure. 

Evidently, due to the existence of maximum for pressure, these black holes enjoy second order phase transition 
in their phase space. The critical pressure and temperature are increasing functions of the massive parameter (Fig. 
m, Cq (Fig. 1111) and dimensions (Fig. 1121) . As for critical horizon radius, it is a decreasing function of the massive 
parameter (left panel of Fig. flOl) and cq (left panel of Fig. fill) while it is an increasing function of dimensions (left 
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FIG. 10: P versus r+ (left panel) and P versus T (right panel) for Co = ci = C 3 = C 4 = 2, C 2 = 3 k = 1, q = 1 and n = 4; m = 0 
(continues line), m — 0.1 (dashed line) and m = 0.2 (dotted line). 




FIG. 11: P versus r+ (left panel) and P versus T (right panel) for ci C3 = C 4 = 2, C2 = 3 k = 1, q = 1, m = 0.2 and n = 4; 
Co = 0.1 (continues line), Co = 1 (dashed line) and Co = 2 (dotted line). 
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FIG. 12: P versus r+ (left panel) and P versus T (right panel) for co = ci = C 3 = C 4 = 2, C 2 = 3 k = 1, q = 1 and m = 0.2; 
n = 4 (continues line), n = 5 (dashed line) and n = 6 (dotted line) in P — r+ diagrams, 
n = 4 (continues line), n = 5 (dotted line) and n = 6 (dashed line) in P — T diagrams. 
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panel of Fig. [13). In addition, the length of coexistence curve is an increasing function of the massive gravity, Cq and 
dimensions. 

The massive parameter, to, is a measurement for the mass of the graviton. Here, we see that by increasing this 
parameter, the critical pressure and temperature increase while critical horizon radius decreases. This indicates that for 
massive graviton, black holes have phase transition in smaller horizon radius. On the other hand, considering massive 
graviton leads to increment of length of coexistence line between two different phases of smaller and larger black holes. 
Therefore, system in this case goes under second order phase transition in higher pressure and temperature. This 
effect highlights one of the differences between gravitational systems with massive gravitons and those with massless 
ones. 


V. CLOSING REMARKS 

In this paper, we have reported a brief discussion regarding to the black hole solutions of massive gravity with a 
linear U (1) gauge field. One of the interesting results of this paper was contribution of the massive gravity to number, 
type and place of horizons. It was shown that the mentioned three properties are highly sensitive to variation of the 
massive parameter and geometrical mass. It was pointed out that considering different values for the parameters will 
lead to different structures for the black holes. Although the massive terms could be up to n + 2 number, it is hard 
to study the total behavior of the solutions with this number of terms. Hence, for the sake of simplicity, we restricted 
our study to the fourth term. 

From thermodynamical point of view, we have studied thermal stability, phase transition and geometrical thermo¬ 
dynamic behavior of the charged black holes in massive gravity. It was shown, in case of finding a critical massive 
parameter for both root and divergence points of the heat capacity, one will come across two sets of equations. In 
case of the root, the presence of all massive coefficients was only observed in the denominator of the obtained relation 
which showed that root of the heat capacity depends on these coefficients. Whereas for the critical massive parameter 
and divergence points of the heat capacity, one of the massive coefficients (ci) had no contribution. 

As for the physical/nonphysical point and second order phase transitions, it was pointed out there are three cases; 
i) only one bound point, ii) one bound point and one second order phase transition, and finally iii) one bound 
point and two second order phase transitions. These three situations are due to the contributions of massive part of 
solutions. In other words, considering the value of massive parameter, one of these cases may happen. Interestingly, 
smaller divergence point of the heat capacity for one dimension coincided with smaller divergence point of the other 
dimensions. The place of coincidences were different for every dimension. In other words, coincidence between two 
different dimensions took place for the specific value of massive parameter. For example, for n = 4 there are different 
values of the massive parameter (say mi, m 2 and etc.) in which the horizon radius for smaller second order phase 
transition will be the same for n = 5, n = 6 and etc. Moreover, upon cancelling the electric charge of the solutions 
this characterized behavior was modified and the presence of second phase transition was vanished. This behavior 
highlights the role of the electric charge in thermodynamical behavior of the system. 

Next, it was shown that thermal stability of the solutions is a function of bound and second order phase tran¬ 
sition points. In other words, considering the number of the phase transition points, system may have different 
stable/unstable states. 

At last, we employed geometrical approach for studying the thermodynamical behavior of the system. It was shown 
that employed thermodynamical metric in this paper, provided a suitable machinery for studying the phase transition 
of these black holes. In other words, regardless of bound point or second order phase transition, all the divergence 
points of Ricci scalar of the considered metric coincided with all root and divergence points of the heat capacity. It 
was also seen that characteristic behavior of the Ricci scalar enables one to recognize the bound point from the second 
order phase transitions. 

In addition, using thermodynamical concept, a study regarding critical behavior of the system was conducted. The 
effects of different parameters on critical values and coexistence curve between two phases were investigated. The 
results regarding the massive parameter were investigated and considerable effects of to on the phase transition status 
were reported. 

The main goal of this paper was investigation of geometrical thermodynamics and van der Waals like phase tran¬ 
sitions in the context of a class of massive gravity. In order to obtain a general model of massive gravity, one may 
follow Stueckelberg mechanism a nd r esuscitate the Goldstone modes which additionally come from breaking diffeo- 
morpliism. For example, in Ref. |l5l| , it was shown that relaxation of the Poincare invariance leads to introduction of 
the large number of the options for generical massive gravity models. It is also worthwhile to consider a dilaton field 
in addition to Maxwell one and study its consequences. Moreover one can apply various nonlinear electrodynamics 
instead of Maxwell field in the context of massive gravity. Generalization of Einstein theory to modified and higher 
derivative gravities with a massive field is another interesting subject. In addition, since obtained black hole solutions 
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have vari ous number and type of horizons, it will be interesting to study thermodynamical relations related to these 
horizons §152j and the possibility of the anti-evaporation property §71, 72]. In addition, following the approach of Ref. 
Ezl , it will be interesting to discuss the effective couplings of action m with scalars to obtain a consistent theory like 
bigravity. All these works could be addressed elsewhere. 
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